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Abstract. We study the evolution of a low-mass X-ray binary by cou- 
pling a binary stellar evolution code with a general relativistic code 
that describes the behaviour of the neutron star. We find that non- 
conservative mass transfer scenarios are required to prevent the formation 
of submillisecond pulsars and/or the collapse to a black hole. We discuss 
the sweeping effects of an active magneto-dipole rotator on the transferred 
matter as a promising mechanism to obtain highly non-conservative evo- 
lutions. 



1. Introduction: a review of the Recycling Scenario 

The "classical" scenario for the formation of millisecond radio pulsar binaries 
with low-mass companions envisages four main stages (see e.g. Bhattacharya & 
van den Heuvel 1991 for a review): (i) The magnetic moment fx of the newly 
formed radio pulsar (/i = BR 3 , where B and R are the surface magnetic field 
and the neutron star, NS, radius 1 , respectively) decays spontaneously on a e- 
folding timescale t^ ~ 10 7 — 10 8 yr (e.g. Lyne, Anderson & Salter 1982) from 
an initial value ~ 10 31 G cm 3 to a final value ~ 10 26 — 10 27 G cm 3 , when the 
decay probably stops (e.g. Bhattacharya & Srinivasan 1986). (ii) Simultaneously, 
the pulsar spin period P increases under magnetic dipole emission according 
to Lpsr = (2/3c 3 )/i 2 (27r/P) 4 (where Lpsr is the bolometric magneto-dipole 
luminosity and c is the speed of light). The pulsar swithces off when eventually 
crosses the "death line" in the \i — P plane, i.e. the line defined by the relation 
^26 /P 2 = 2 x 10 3 (where /^26 = ^/10 26 G cm 3 ) below which it is believed that the 
radio pulsar phenomenon does not take place (see e.g. Ruderman & Sutherland 
1975). (iii) The companion star overflows its Roche lobe and transfers mass with 
angular momentum to the NS via a Keplerian accretion disc, thereby spinning 
it up to millisecond periods (close to the Keplerian period at inner rim of the 
accretion disk, see below) and back across the death line (recycling). During 



1 Actually, since the NS is a relativistic object, R is the NS circumferential equatorial radius, 
i.e. the proper circumference in the equatorial plane divided by 2n. 
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this phase the system is visible as a low-mass X-ray binary (LMXB). (iv) Mass 
transfer ceases: the NS is again visible as a radio pulsar whose spin rate decays 
under magnetic dipole emission very slowly as Lpsr oc /i 2 and fi is reduced 
by three or four orders of magnitude. The end point is therefore a millisecond 
pulsar orbiting a low mass companion (< 0.3 M ) that is the remnant of the 
~ 1 M Q mass donor. 

Recycling in Transient Systems 

Most LMXBs are NS Soft X-ray Transients (NSXT), i.e. transient systems har- 
boring a NS (see Campana et al. 1998 for a review). Adopting the same conver- 
sion efficiency of the accreting matter energy into X-rays during the outbursts 
and the quiescent states (but see Barret et al. 2000 for a different explanation), 
the inferred variations in the accretion rate £1X6 cl factor ~ 10 . 

NSXTs can provide a direct evidence of the recycling scenario, since during 
stage (iii) the mass transfer rate varies up to five orders of magnitude. During 
the LMXB phase the accretion disk is truncated because of one of the following 
reasons: (i) the interaction with the magnetic field of the NS, which truncates the 
disc at the magnetospheric radius Rm, at which the accretion flow is channeled 
along the magnetic field lines towards the magnetic poles onto the NS surface; 
(ii) the presence of the NS surface itself at R; and (iii) the lack of closed Keplerian 
orbits for radii smaller than the marginally stable orbit radius, i?MSO ( a t few 
- depending on the mass and spin of the compact object - Schwarzschild radii 
from the NS centre). The position of Rm is determined by the istantaneous 
balance of the pressure exerted by the accretion disc and the pressure exerted 
by the NS magnetic field: 

Rm = 1.0 x 10 6 <p fi% 7 m- 1/7 RQ 2/7 m- 2/r cm , (1) 

where <f> < 1, m is the NS gravitational mass in M 2 , R§ is the NS radius in 
units of 10 6 cm, and rh is the baryonic mass accretion rate 2 in Eddington units 
(the Eddington accretion rate is 1.5 x 10 _8 i?6 M Q yr _1 ). Eq. 1 shows that as 
m decreases, Rm expands. 

Accretion onto a spinning magnetized NS is centrifugally inhibited once 
Rm expands beyond the corotation radius i?co> & t which the Keplerian angu- 
lar frequency of the orbiting matter is equal to the NS spin: Rqo = 1-5 x 
10 m 

l/3p2/3 

cm where P_3 is the NS spin period in milliseconds. In this case 
the accreting matter could in principle be ejected from the system: this is called 
propeller phase (Illarionov & Sunyaev 1975). Finally, if Rm further expands be- 
yond the light-cylinder radius (where an object corotating with the NS attains 
the speed of light, i?LC = 4.8 x 10 6 P_3 cm), the NS becomes generator of mag- 
netodipole radiation and relativistic particles. Indeed, a common requirement 
of all the models of the emission mechanism from a rotating magnetic dipole is 
that the space surrounding the NS is free of matter up to .Rlc- 



2 Actually, since the NS is a relativistic object, it is important to distinguish between the baryonic 
mass, roughly speaking a measure of the amount of matter, and the gravitational mass that is 
smaller by a factor of ~ (3/5)GM 2 / R ~ O.lMc 2 , corresponding to the binding energy of the 
NS. 
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Let us consider the behaviour of a NSXT at the end of an outburst. Adopt- 
ing m ~ 1 in outburst, eq. 1 gives Rm ~ R ~ 10 6 cm. In quiescence m ~ 10~ 5 
and Rm = (1(T 5 )~ 2 / 7 x 10 6 = 2.7 x 10 7 cm > .Rlc, for spin periods up to few 
milliseconds. Therefore it is likely that, during the quiescent phase, a magneto- 
dipole emitter switches on (see e.g. Stella et al. 1994; Burderi et al. 2001). In 
this the NS fj, and P place such a system above the "death line", it is 

plausible to expect that the NS turns-on as a millisecond radio pulsar until a 
new outburst episode pushes Rm back, close to the NS surface, quenching radio 
emission and restoring accretion. 

2. To accrete: conservative mass transfer 

A rapidly rotating NS can support a maximum mass (against gravitational col- 
lapse) much higher than the non-rotating mass limit, since the centrifugal force 
attenuates the effects of the gravitational pull. Conversely, if a rotating NS has 
a mass that exceeds the non-rotating limit (supramassive NS), it will be subject 
to gravitational collapse if it loses enough angular momentum, J. In contrast to 
the standard behaviour, supramassive NSs spin-up just before collapse, even if 
they lose energy (Cook, Shapiro & Teukolsky 1992). The value of the maximum 
rotating and non-rotating mass depends on the equation of state (EOS) govern- 
ing the NS matter, and the minimum allowed period for a given mass and EOS 
occurs when gravity is balanced by centrifugal forces at the NS equator (mass 
shedding limit). In this case the NS has the maximum J allowed for that given 
mass, J max - The NS radius, which is always in the order of 10 6 cm, depends on 
the mass, EOS, and J of the NS: a rapidly rotating NS can have a much larger 
radius than a non-rotating one (up to 40%, Cook, Shapiro k, Teukolsky 1994). 

In order to study the evolution of a LMXB in detail we coupled a modi- 
fied version of rns (rotating NSs) public domain code by Stergioulas & Morsink 
(1999), that allowed us to build a complete grid of relativistic equilibrium config- 
urations, with the binary evolution evolution code ATON (D'Antona, Mazzitelli 
k Ritter, 1989). We assume the NS to be low-magnetized (B ~ 10 8 G). Our 
computations show that during the binary evolution, the companion transfers as 
much as 1 M Q to the NS, with an accretion rate of ~ 10~ 9 M yr _1 . This rate is 
sufficient to keep the inner radius of the accretion disc, rjj, in contact with the 
NS surface, thus preventing the onset of a propeller phase capable of ejecting a 
significant fraction of the matter transferred by the companion. 

Matter leaves the inner radius of the disk td with a specific angular mo- 
mentum j that is a function of ro, Mq, and J (Bardeen 1970) We can therefore 
write the evolutionary equations: 

dM B /dt = M c 
dJ/dt = f(J) ■ M B 

where Mb is the baryonic mass of the star and Mc is the mass transfer rate 
from the companion, and /(J) = j for J < Jmax (Mb) and /(J) = dJ max /dM B 
for J = Jmax (Mb), since we have assumed that at mass shedding the matter of 
the disc dissipates angular momentum in order to accrete onto the NS, keeping 
the NS at the mass shedding limit. 



(2) 
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Figure 1. Left: Evolution of the NS in the mass-radius plane. La- 
belled dots indicate the corresponding value of the NS spin period in 
ms. Right: Radial dependence of the pressures relevant for the evolu- 
tion of accreting NSs and recycled pulsars. The parameters adopted 
are //26 = 5, P-;$ = 1.5, Rq = 1, m = 1.4 

Once the accretion ends, Rm expands beyond i?LC and a magneto-dipole 
rotator (radio pulsar) switches-on. In this situation is convenient to describe the 
NS evolution in terms of baryonic mass and total energy rather than in terms 
of baryonic mass and angular momentum: 



We integrate the differential equations 2 or 3 using a finite-differences method 
(see Lavagetto et al. 2004 for details). 

In Figure 1 (left panel) we plot, in the mass-radius plane, the evolution of 
a system with a donor star of 1.15M0, which loses 0.91M© during the accretion 
phase, which is representative of most of the simulations. To describe the NS 
we adopted the quite standard EOS FPS. The upper dashed curve is the mass- 
shedding limit, the lower dashed curve is the the sequence of stable non-rotating 
configurations, the thin solid line on the lower right is the limit beyond which 
the NS collapses onto a black hole, the thick solid line shows the evolution of 
the NS during the accretion phase, and, finally, the thick dashed line shows the 
subsequent radio-pulsar phase. The first phase of accretion is charachterized by 
a rapid spin up that brings the star to the mass-shedding limit, along which it 
evolves with P slowly decreasing from 0.78 to 0.58 ms. Once the accretion stops 
during the radio-pulsar phase, P decreases (because the NS is supramassive) to 
0.54 ms, until the NS eventually collapses to a black hole. 

All of our simultations indicate that: a) for most EOS an accretion induced 
collapse to a black hole is almost unavoidable; b) in the absence of an external 
mechanism that brakes down the NS (such as gravitational waves), a small 
fraction of the accreted matter (~ O.1M ) is sufficient to spin-up the NS to 
periods below one millisecond. In contrast, no submillisecond pulsars have been 



f dM B /dt 
\ dM G /dt 



= 



= -(2/3c 5 )^ 2 (27r/P) 4 
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detected up to date: the shortest observed spin period is 1.5 ms (Backer et 
al. 1982), uncomfortably higher than the theoretical predictions. Both these 
predictions are in contrast with the existence of the population of MSP all 
spinning at frequencies well above 1 ms. 

3. Not to Accrete: non conservative mass transfer 

A non-conservative mass transfer has been often invoked to overcome the diffi- 
culties outlined above. Indeed the ejection of more than 90% of the transferred 
mass could prevent the gravitational collapse and keep the NS spin period above 
1 ms. The propeller effect discussed in §1 is, in principle, a promising mechanism 
to obtain highly non-conservative evolutions. However, the virial theorem sets 
stringent limits on the fraction of matter that can be ejected in this phase. In 
fact, it states that, at any radius in the disk, the virialized matter has already 
liberated (via electromagnetic radiation) half of its available energy. Considering 
that ro ~ R for these kind of systems, the accreting matter has already radiated 
50% of the whole specific energy obtainable from accretion (e acc ~ GMq/R) be- 
fore setting onto the NS surface. Therefore to eject this matter ~ 0.5e acc must 
be given back to it. As the only source of energy is that stored in the NS as 
rotational energy by the accretion process itself, the typical ejection efficiency 
is < 50%, considering that once the system has reached the spin equilibrium, 
no further spin-up takes place and the storage of accretion energy in rotational 
energy is impossible. Thus, the accreted mass is > 0.5 M Q , well enough to 
spin-up the NS to submillisecond periods. An alternative viable hypothesis to 
explain the lack of ultrafast rotating NSs is that gravitational wave emission 
balances the torque due to accretion (see a review in Ushomirsky, Bildsten, & 
Cutler 2000), although this effect cannot prevent the gravitational collapse. 

The only way to overcome these difficulties is to obtain ejection efficiencies 
close to unity. This is indeed possible if the matter is ejected so far away from 
the NS surface (r >> R) that it has an almost negligible binding energy e = 
GMc/r « e acc . As the NS is spinning very fast, the switch-on of a radio 
pulsar is unavoidable once Rm > Rlc- in this case the pressure exerted by the 
radiation field of the radio pulsar may overcome the pressure of the accretion 
disk, thus determining the ejection of matter from the system. Once the disk has 
been swept away, the radiation pressure stops the infalling matter as it overflows 
the inner Lagrangian point L±, where e ~ 0. 

The push on the accretion flow exerted by the magnetic field of the NS can 
be described in terms of an outward pressure (we use the expressions outward or 
inward pressures to indicate the direction of the force with respect to the radial 
direction): Pmag = B 2 /4:TT = 7.96x 10 14 ^26 r 6~ 6 dy/cm 2 , where r$ is the distance 
from the NS center in units of 10 6 cm. If the disk terminates outside Rlc, 
the outward pressure is the radiation pressure of the rotating magnetic dipole, 
which, assuming isotropic emission, is Ppsr = 2.04 x 10 12 PZ% ^26 r 6~ 2 dy/cm 2 . 
In Figure 1 the two outward pressures (Pmag and -Ppsr) are shown as bold lines 
for typical values of the parameters (see figure caption). 

The flow, in turn, exerts an inward pressure on the field. For a Shakura- 
Sunyaev accretion disc (see e.g. Frank, King & Raine 1992): Pqisk = 1-02 x 

10 l 6Q ,-9/10 n -l i5L 17/20 m i /4 ^17/20^17/5^21/8 dy/cu] 2^ where Q ig thg Shakura _ 
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Sunyaev viscosity parameter, no.615 = n/0.615 ~ 1 for a gas with solar abun- 
dances (where n is the mean particle mass in units of the proton mass), and 
/ = [1 - {Re/re) 1 / 2 } 1 / 4 < 1. We measure M in units of L 37 = L/10 37 ergs/s 
from L = GMM/R. 

In Figure 1 (right panel) the inward disc pressure for a luminosity Lmax> 
corresponding to the outburst luminosity of an NSXT, is shown as a thin solid 
line. The disc pressure line, which intersects -Pmag at Rlc, defines a critical 
luminosity Switch) shown as a dashed line in Fig. 1, at which the radio pulsar 
switches-on. 

The intersections of the Pdisc line corresponding to Lmax with each of 
the outward pressure lines define equilibrium points between the inward and 
outward pressures. The equilibrium is stable at r = R m , and unstable at 
r = i?sTOP ; which can be derived equating Ppsr and -Pdisk: -Rstop ~ 8 x 

10 11 «^6/ 2 5 no J/^34/ 2 5 /136/25L 34/ 2 5 mV25 ^ 26 16/5 p 3 23 /5 ^ ^ ^ &g ^ . g 

steeper than Pqisc j if a small fluctuation forces the inner rim of the disc inward 
(outward), in a region where the magnetic pressure is greater (smaller) than the 
disc pressure, this results in a net force that pushes the disc back to its original 
location R m . As Ppsr is flatter than -Pdisc> with the same argument is easy to 
see that no stable equilibrium is possible at -Rstop and the disc is swept away 
by the radiation pressure. This means that, for r > i?sTOP> no disc can exist for 
any luminosity < Lmax- The sudden drop in the mass-transfer rate during the 
quiescent phase of a NSXT initiates a phase that we termed "radio ejection" , in 
which the mechanism that drives mass overflow through L\ is still active, while 
the pulsar radiation pressure prevents mass accretion. As the matter released 
from the companion cannot accrete, it is now ejected as soon as it enters the 
Roche lobe of the primary. The distance of L\ from the NS center, i?Li, depends 
only on the orbital parameters. In the approximation given by Paczynski (1971), 
we can impose -Rstop/-Rli = 1 and solve for the orbital period: 

P crit = 1.05 x (a- 36 n ^? 5 i?r) 3 / 50 L^/ 25 rnVio^ 26 24/5 p 48/5 5(m) ^ fa (4) 

where g(m, 1712) = [1 — 0.462 (7712 /(m + ttt^)) 1 / 3 ] -3 / 2 ^ + m^) -1 / 2 , and 777-2 is the 
mass of the companion in solar masses. For a system with P or b > -Pcrit (e.g. 
-forb = 3 h in Fig. 1), once a drop of the mass-transfer rate has started the radio 
ejection, a subsequent restoration of the original mass-transfer rate is unable to 
quench the ejection process (see Burderi et al. 2001 for a detailed discussion of 
this effect), thus prolonging the ejection phase. 

This "radio-ejection" phase envisaged by Burderi et al. (2001), has been 
spectacularly demonstrated by the discovery of PSR J1740-5340, an eclipsing 
MSP, with a spin period of 3.65 ms and an orbital period of 32.5 h (D'Amico 
et al. 2001), very close to P cr it for the parameters of this system (Burderi, 
D'Antona, & Burgay 2002). The peculiarity of this system is that the companion 
is still overflowing its Roche lobe. This is demonstrated by the presence of 
matter around the system that causes the long lasting and sometimes irregular 
radio eclipses, and by the shape of the optical light curve, which is well modeled 
assuming a Roche-lobe deformation of the mass losing component (Ferraro et al. 
2001). An evolutionary scenario for this system has been proposed by Burderi 
et al. (2002), who provided convincing evidence that PSR J1740-5340 is an 
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example of a system in the radio-ejection phase, by modeling the evolution of 
the possible binary system progenitor. 

In conclusion, while the evolution without a radio-ejection phase implies 
that a large fraction of the transferred mass is accreted onto the NS (because of 
the constraints imposed by the virial theorem), we have demonstrated that the 
switch-on of a radio pulsar (associated to a significant drop in mass transfer) 
could determine long episodes characterized by ejection efficiencies close to 100%, 
as the matter is ejected before it falls into the deep gravitational potential well of 
the primary. This phenomenon can prevent gravitational collapse and ultra-fast 
spinning NSs. 
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